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ABSTRACT 

Context. Diffuse interstellar bands (DIBs) are non-stellar weak absorption features of unknown origin found in the spectra of stars 
viewed through one or several clouds of the interstellar medium (ISM). Research of DIBs outside the Milky Way is currently very 
limited. In particular, spatially resolved investigations of DIBs outside of the Local Group are, to our knowledge, inexistent. 

Aims. In this contribution, we explore the capability of the high-sensitivity Integral Field Spectrograph, MUSE, as a tool for mapping 
diffuse interstellar bands at distances larger than 100 Mpc. 

Methods. We used MUSE commissioning data for AM 1353-272 B, the member with the highest extinction of the Dentist’s Chair, an 
interacting system of two spiral galaxies. High signal-to-noise spectra were created by co-adding the signal of many spatial elements 
distributed in a geometry of concentric elliptical half-rings. 

Results. We derived decreasing radial profiles for the equivalent width of the 45780.5 DIB both in the receding and approaching side 
of the companion galaxy up to distances of ~4.6 kpc from the centre of the galaxy. The interstellar extinction as derived from the 
Hff/H/3 line ratio displays a similar trend, with decreasing values towards the external parts. This translates into an intrinsic correlation 
between the strength of the DIB and the extinction within AM 1353-272 B, consistent with the currently existing global trend between 
these quantities when using measurements for Galactic and extragalactic sightlines. 

Conclusions. It seems feasible to map the DIB strength in the Local Universe, which has up to now only been performed for the 
Milky Way. This offers a new approach to studying the relationship between DIBs and other characteristics and species of the ISM in 
addition to using galaxies in the Local Group or sightlines towards very bright targets outside the Local Group. 
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1. Introduction 


Diffuse interstellar bands (DIBs) are non-stellar weak absorp¬ 
tion features found in the spectra of stars viewed trough one or 
several clouds of the interstellar mediumjTSM, Herb igj 119951) . 
They were identified for the first time bv lHeger , ( 1922jh but their 
interstellar origin was established in the 1930s (Me rrilliri934l) . 
Almost one century after their discovery, the nature of their car- 
riers (i.e. the agent that cau ses these features) remains a mystery 
(iFulara & Krelowskil 120001) . However, several of the DIBs are 
relatively strong and present good correlations with the amount 
of neutral hydrogen, the extinction, and the interstellar Na I D 
and Ca H and K lines alon g a given line of sight (e.g. iHerbigl 
1 19931: [Frie dman et al. 2011). Thus, irrespective of the actual na¬ 
ture of carriers, DIBs can be used as tools to infer the prop¬ 
erties of the ISM struct ure, as shown by several recent works 
for the Milky Way (e.g. Pusgitarinietakl 20151 IZasowsk i et al.1 
l20lllKoset al.ll2014l:ICordiner et al.ll2013l). 


* Based on observations made with ESO telescopes at the La Silla 
Paranal Observatory under program ID 60.A-9100(B) 


Equivalent examples for other galaxies are difficult to find. 
There is a certain level of spatial resolution in the inve stigatio ns 
of DIBs in galaxies within the Local Group (|Ehrenfreund et alJ 
l2002t ICordiner et akll201 It Ivan Loon et al.H20lT 7 since several 
sightlines in a given galaxy are sampled. Works targeting DIBs 
outside the Local Group are still rare and address individual 
sightlines for a gi ven target (e.g. iHeckman & Lehnertl [2000; 
[Phillips et al.ll2013l) . Currently, no mapping of a DIB has been 
performed for galaxies outside of the Local Group (beyond ~ 1 
Mpc). 

The advent of highly efficient integral field units renders this 
feasible for the first time. They simultaneously record the spa¬ 
tially resolved spectral information in an extended continuous 
field. This can be tessellated a posteriori, preserving the rele¬ 
vant spatial information as best possible, but at the same time 
optimizing the quality of the signal by co-adding the spectra be- 
longing to a given portion of the mapped ar ea - the tile - (e.g. 
iMonreal-Ibero et aTl2005l:lSandin et al.ll2008l) . Here, we use this 
principle to explore the capabilities of MUSE to obtain spatially 
resolved information of extragalactic DIBs. 
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We carried out our experiment in the interacting system 
AM 1353-272. The system consists of two components. The 
main galaxy (A) presents two promi nent ~40 kpc long tidal 
tails and w as thoroughly studied by Weilbac her et alJ (2000, 
120021 l2003h . The companion (B) is a low-luminosity (Mg = 
-18.1 mag) disk-like galaxy of disturbed morphology undergo¬ 
ing a strong starburst and with high extinction. The relati ve posi - 
tions of A and B are shown in Fig. 1 of lWeilbacher et akl (120021) . 
Here, we used AM 1353-272 B as a test bench to explore the 
possibility of DIB detection and mapping. We assume a distance 
of D = 159 Mpc for AM 1353-272 using Hq = 75 km s _I , as in 
[Weilbacher et al.l (i2002l) . This implies a linear scale of a 0.771 
kpc arcsec -1 . 


2. Observations and data reduction 

The interacting system AM 1353-272 was observed as part of 
commissioning run 2a dBacon et alJl201 4l) of the MUSE instru¬ 
ment at the VLT. Six 900 s exposures were taken in two posi¬ 
tions between 2014-04-29T04:22:00 and 06:00 under conditions 
of a few thin clouds and an average effective airmass of 1.04. 
The seeing as given by the Gaussian FWHM of the stars in the 
white-light image of the final cube is 078. Each position was 
observed at position angles (P.A.) of 0°, 90°, and 180°. The in¬ 
strument was set to wide-held mode (~ 1 ' x 1 ' held of view) with 
a nominal wavelength range (~ 4800.. .9300A). The instrument 
resolution at 7000A is R ~ 2800. The standard star GJ 754.1 A 
was observed on 2014-04-28T09:57:00, and the reduction made 
use of twilight skyhats taken on 2014-04-27. 

Standard reduction steps were followed with MUSE pipeline 
vl.O ( [Weilbacher et al.l2014l) . with bias subtraction, hat-fielding, 
and spectral tracing using the lamp-hat exposures, wavelength 
calibration, all using standard daytime calibration exposures 
made the morning after the observations. The standard geometry 
table and astrometric solution for commissioning run 2a were 
used to then transform the science data into pixel tables. This 
last step included correction for atmospheric refraction, flux cal¬ 
ibration, and sky subtraction. The sky spectrum was measured in 
the blank areas of each exposure, set to the darkest 40% of the 
held of view, and data were corrected to barycentric velocities, 
resulting in corrections of -0.6 to -0.7 km s _I for the time of ob¬ 
servations. The hnal, combined datacube was then reconstructed 
from all six exposures, sampling the data at 072 x 072 x 0.8 A, 
rejecting cosmic rays in the process. We used a subset of 111x91 
spatial elements for the experiment presented here. 


3. Extracting information 

3. 1. Tile definition 

The galaxy has receding and approaching velocities at the south¬ 
east and northwest side, respectively (Weilbacher et al. in prep.). 
Since the velocity gradient is strong (Ai’>300 km s 1 ), we treated 
each side of the galaxy independently, this way avoiding a large 
artificial broadening of the spectral features. To do this, the tiles 
were defined as elliptical half rings following the morphology of 
the galaxy, with a fixed P.A.= 153°, measured counterclockwise, 
relative to the north, and an ellipticity of b/a = 0.2. The ring 
sizes represent a compromise between retaining as much spatial 
resolution as possible and having a reliable DIB detection. The 
final tiles (depicted in Fig. [T]) have between 29 and 159 spatial 
elements each (projected areas between 0.7 and 3.8 kpc 2 ). 


Fig. 1. Apertures used to 
extract the high signal- 
to-noise spectra (i.e. the 
tiles). Each tile has been 
coloured differently using 
a palette that follows the 
blue-to-red velocity distri¬ 
bution within the galaxy. 
As a reference, the re¬ 
constructed white-light im¬ 
age is overplotted as con¬ 
tours in logarithmic stretch¬ 
ing with 0.25 dex steps. 


3.2. Feature fitting 

The spectra belonging to each tile were co-added, extracted, 
and corrected from a Galactic extinction of A v = 0.165 
dSchlaflv & Finkbeinerl f20Tl ) using the IRAF _task deredden 
and a ssuming Ry = A(V)/E(B — V) = 3.1 dRieke & Febofskvl 
119851) . In each spectrum, we fitted different spectral features 
needed for our analysis using the IDE-based routine mpfit 
(lMarkwardti:2()09 ). The results discussed here are based on the 
equivalent width of the DIB at 45780.5 (EW(45780.5)), and 
fluxes for Ha, H8, rOnf|4 5007, and [Nn]46584. We refer to 
[Monreal-Ibero et ahl ( 2010 ) for a description of the information 
recovery for the emission lines since it is quite standard and not 
particularly challenging for the experiment described here. 

The most critical feature is the DIB at 45780.5. This was fit¬ 
ted together with the sodium doublet at 445889.9,5895.9 and the 
Hei45876 emission line using Gaussian functions, fixing their 
relative central positions according to a redshift of z = 0.039467 
(v = 11 840 km s _1 ), as determined from the average of the 
Ha centroid in the two inner half rings and in accord with the 
11 790 + 50 kms~' reported bv lWeilbacher et aD (120031) and with 
same line width. Additionally, we included a one-degree polyno¬ 
mial to model the stellar continuum. The DIB at 45797.1 was not 
included in the fit since we did not see any evidence for it in our 
spectra. It is known that the ratio of the equivalent width of these 
two DIBs can vary up to a factor ~5, depending on the charac¬ 
teristics of the ISM in the sight line of consideration (lYos et al.l 
[201 lh . but only in very specific situations both DIBs are compa¬ 
rable. Typically, the DIB 45797.1 is ~3-5 times fainter. 

Figure [2] shows the six fits for this critical range. The DIB 
at 45780.5 is clearly seen in all of them, with the possible ex¬ 
ception of the outermost approaching half-ring. Hereafter, error 
bars indicate a prelimin ary estimate of the uncer tainties as de¬ 
rived using Eq. 1 in iMonreal-Ibero et ah ! (2 01 lb . which takes 
into account the signal-to-noise ratio of the spectra. However, 
the greatest source of uncertainty is the possible contamination 
in the spectra by a stellar spectral feature at 45782, refff frame, 
with contributions from Fei, Cri, Cui, and Mgi dWorthev et al.l 
119941) . and whose strength increases with age and metallicity. 
An in-depth study to model this feature is beyond the scope of 
this initial experiment. In the following we discuss the expected 
maximum amount of contamination and provide an ad hoc cor¬ 
rection for it. 

We measured EW(Ha) in emission varying from ~ 90 A in 
the inner rings to ~ 60 A in the outer ones, which points to 
a very young (i.e. <7 Myr) underlying stellar population. The 
ages of the star-formin g knots in AM 1353-272 A (<40 Myr, 
IWeilbacher et al.ll2000l) , as indicative of the moment when the 
interaction triggered the star formation would be another way 
to estimate the age of the stars in AM 1353-272 B. Fikewise, 
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Fig. 2. Spectra for the six half rings showing the fit in the area of the 
DIB, ordered from smallest top to largest bottom velocity and normal¬ 
ized to the median in the displayed spectral range. The expected loca¬ 
tions of the discussed features are marked at the bottom of the figure. 
The colour code in the fits echoes the ordering in velocity. Likewise, the 
shift of the spectral features is clearly visible. 


using the N 2 and (93A2_calibrators and the expressions pro¬ 
vided by [Marino et al.1 (12013l) . we estimate a metallicity for the 
inner rings of 12 + log(<9///) = 8.50 and 8.42, respectively 
(i.e. about half solar), decreasing outwards. To estimate an up¬ 
per limit for the contamination of the stellar feature, we con¬ 
servatively created 1000 mock spectra by taking the spectrum 
of a 60 Myr single s tellar population from the M ILES library 
with solar metallicity dFalcon-Barroso et al.l201 lb . adding Gaus- 
sians representing the DIBs at /15780.5 and T5797.1 and the 
He iT5876 emission line. Then, we added noise with a scale of 
the standard deviation between 5 020 A and 5 040 A, rest-frame. 
EWL15780.5) in the spectra was varied between 10 mA and 
735 mA. As discussed above, a common range for the rela¬ 
tive strength of the two DIBs is ~ 3-5. Here, we assumed 
EW(/15797.1) = EW(/15780.5)/3. We measured EW(T5780.5) in 
these mock spectra using the same methodology as in the real 
data. The relation between input and output equivalent widths 
was fitted to a one-degree polynomial that was used to esti¬ 
mate an upper limit to the correction due to the stellar fea¬ 
ture contamination. Measured EW(715780.5)<200mA are con¬ 
sistent with no DIB, EW(T5780.5)~ 300 mA should be reduced 
by a 40%. For EW(/l5780.5)>500niA a reduction of ~20% is 
enough. This simulation also suggests that under typical con¬ 
ditions, an EW(/t5780.5)>800 mA would be necessary to have 
some expectations for detecting the DIB at T5797.1. 


4. Results 

4.1. Radial profiles 

The upper plot of Fig. [3] shows EW(715780.5) as a function of 
the semimajor axis. The receding and approaching halves of the 
galaxy both present a negative gradient of the DIB strength with 
distance to the centre. This points towards a direct link with the 
intrinsic properties of the galaxy. 


One of the best established correlations between DIB 
strength and other characteristics of the ISM is that with the 
interstellar extinction dHerbid fl993h . We derived an e stima te 
of the extinction using the extinction curve of iFluks et all 
(fl994t) and assuming an intrinsic Balmer emission-line ratio 
H a/ Hj3 = 2 .86 for a case B aproximation and T e = 10000 K 
(fOsterbrock & Ferlandll2006h . The corresponding radial profiles 
are shown in the lower plot of Fig. [3] As with the DIB, they dis¬ 
play a decreasing gradient towards larger distances from the cen¬ 
tre of the galaxy. This implies an intrinsic correlation between 
extinction and EW(/15780.5) for AM1353-272 B and further 
supports the derived radial profiles for the DIB. 



Fig. 3. Radial profiles for EWL15780.5) {top) and the extinction ( bot¬ 
tom) as derived from Ho'/H/f for the approaching (stars and blue lines) 
and receding (circles and red lines ) sides. The horizontal bars traversing 
the filled symbols mark the radial extent over which the spectra were co¬ 
added. In the upper plot, filled symbols indicate measurements without 
correction, while open symbols include the maximum correction for the 
contamination of the stellar feature at T5782 A. 


4.2. Comparison with other sightlines 

Figure [T] illustrates the comparison between this intrinsic cor¬ 
relation with measurements for other sightlines. The sample is 
by no means complete, but covers a diversity of environments, 
methodologies, sampled areas, etc. This includes examples of 
spatially resol ved measurements of the two large spirals i n the 
Local Group dPuspitarini et al.1 [2 0131: ICordiner et ~akl 1201 lit as 
well as for the Magellanic Clou ds (IW eltv et ak|2006l). Likewise, 
we inclu ded some supernovae ( Phillips et al.1 2013 ) and nearby 
galaxies dHeckman & Lehnertl l2000l) . There is a clear correla- 
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tion between extinction and EW(/15780.5). Specifically, the ratio 
between this two quantities across AM 1353-272 B is constant 
(within experimental errors), despite variations in metallicity, ra¬ 
diation field, etc. (Weilbacher et al. in prep). When examining 
the ensemble of data globally, the dispersion is larger, pointing at 
secondary parameters modulating this correlation. In particular, 
for a given reddening, DIBs in the Magellanic Clouds are ajactor 
2-6 fainter than Galactic sightlines, as pointed out bv lWeltv et al.i 
(120061) . All of our measurements in AM 1353-272 B agree well 
with this general relation. Interestingly enough, the relation be¬ 
tween EW(T5780.5) and E(B - V ) is tighter when only resolved 
studies of disk-like galaxies are taken into account (Milky Way, 
Andromeda, our data). AM 1353-272 B indeed occupies an ex¬ 
tension towards high EWs and E(B-V) of the relation found for 
the large spirals in the Local Group. 



l°g(E(B—V)) 

Fig. 4. Relation between EW(d5780.5) and extinction for AM1353-272 
(red circles and blue stars). Symbols and colours are as in Fig. [3] Ad¬ 
ditionally, we show the following Galactic (triangles ) and extragalac- 
tic (in verted triangl es) sightlines: green (Milky Way, Pusnitarini et all 
20131); purple (M 31, ICordiner et al. fed 111) . orange (Magel lanic Clouds, 
Weltv et all 12006T), light blue f supernovae. | Phillins et alll20l 3l) . black 
(dusty starbursts, Heckman & Lehnert 200CJ). We used filled symbols 
when a galaxy had spatially resolved detections and open symbols when 
there is a single sightline for a given target. 


5. Summary and perpectives 

This work presents the first spatially resolved detection of a DIB 
in a galaxy outside the Local Group. By using 568 MUSE spec¬ 
tra sampling a total projected area of 13.5 kpc 2 , we were able 
to measure the equivalent width of the DIB at /15780.5 in six lo¬ 
cations of AM 1353-272B. Our strategy constitutes an alterna¬ 
tive approach to extragalactic DIB studies different from using 
targets in the Local Group and individual sightlines outside the 
Local Group. We found three main results: 

1. We found decreasing radial profiles for EW(T5780.5), both 
in the receding and approaching sides of AM 1353-272 B up 
to distances of ~4.6 kpc from the galaxy centre. 

2. The interstellar extinction displays a similar trend, with de¬ 
creasing values towards the external parts. This translates 
into a correlation between the strength of the DIB and the 
extinction within AM 1353-272 B. 

3. A comparison of E(B - V ) and EW(/15780.5) in AM 1353- 
272 B and other sightlines shows that this intrinsic corre¬ 
lation agrees with the existing global trend between these 
quantities, especially when compared with targets in the spi¬ 
ral galaxies of the Local Group (i.e. our Galaxy and M31). 


When seen globally, these three results demonstrate that the 
spatially resolved detection of DIBs in galaxies outside the Lo¬ 
cal Group is possible and 2D DIB mapping feasible thanks to 
high-sensitivity instruments like MUSE. There are currently no 
2D maps of DIBs in galaxies other than the Milky Way. The 
construction of such maps will permit studying the relationship 
between the DIBs and other components of the ISM in condi¬ 
tions different from those found in our Galaxy. We will be able 
to address questions like whether an extinction-DIB correlation 
in extremely dense environments also exists, to which degree the 
strength of the different DIBs depends on the dust-to-gas ratio, 
how DIBs react to the radiation field typical of an environment of 
extreme star formation, or whether the different relations locally 
found (e.g. the various DIB families) also exist in other galax¬ 
ies, and how they depend on the intrinsic characteristics of the 
galaxies (e.g. morphology, metallicity, etc.). The answers will be 
a critical test that will help to constrain the nature of the carriers. 
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